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The extraordinary properties of resonant four-wave mixing of backward waves in doped negative- 
index materials are investigated. The feasibility of independent engineering of negative refractive 
index and nonlinear optical response as well as quantum control of the nonlinear propagation process 
in such composites is shown due to the coherent energy transfer from the control to the signal 
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amplification, and realizing a miniature mirrorless optical parametric generator of the entangled 
backward and ordinary waves are among the possible applications of the investigated processes. 
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INTRODUCTION 

Optical negative index metamaterials (NIMs) with si- 
multaneously negative electric and magnetic responses 
(NIMs) promise revolutionary breakthroughs in pho- 
tonics (see, for example, [1]). The most detrimen- 
tal obstacle toward the application of NIMs is strong 
the absorption of light that is inherent to this class of 
materials. The possibility to overcome such obstacles 
based on three- wave [x^j optical parametric amplifi- 
cation (OPA) in NIMs was proposed in 0, 0). How- 
ever, a strong nonlinear-optical (NLO) response must be 
ensured and the corresponding frequency domain must 
overlap the phase-matching frequency-interval to realize 
such feasibility. Nonlinear optics in NIMs remains so 
far a less-developed branch of optics. The possibility of 
quadratic NLO optical response in NIMs attributed to 
the asymmetry of the voltage-current characteristics of 
their nanoscaled building blocks was predicted in 0, Q • 
Recent experimental demonstration of the exciting feasi- 
bilities to craft NIMs with strong NLO response has been 
reported in Q. Counterintuitive properties of nonlinear 
propagation processes in NIMs with NLO response, 
such as second harmonic generation and three- wave OPA, 
as compared with their counterparts in natural materials, 
were revealed in @, i, 0, d, H . The striking changes in the 
optical bistability in layered structures including MIM 
were shown in 10]. A review of the corresponding theo- 
retical approaches is given in . Extraordinary proper- 
ties of a three- wave mixing (TWM) backward-wave (BW) 
optical parametric oscillator (OPO) were predicted sev- 
eral decades ago in HjlljHl- However, phase matching 
of the coupled waves with the opposite orientation of their 
wave vectors, which is required for mirrorless oscillations 
in ordinary materials, is very difficult to achieve for the 
waves with substantially different frequencies. For the 
first time, TWM backward-wave mirrorless optical para- 
metrical oscillator (BWMOPO) with all three different 
optical wavelengths was realized only recently [HI, [l6| . 
Phase matching of waves with an antiparallel orienta- 



tion of their wave vectors has been achieved in a submi- 
crometer periodically poled NLO crystal, which has be- 
come possible owing to recent advances in nanotechnol- 
ogy. As outlined, a major technical problem in creating 
BWMOPO stems from the requirement of phase match- 
ing for traveling waves with oppositely oriented wave vec- 
tors, a situation that is inherent to ordinary materials. 
However, the situation dramatically changes in NIMs. 
Due to the opposite orientation of the energy flow and 
wave vector, which is inherent to NIMs, phase match- 
ing of backward waves becomes possible for the parallel 
orientation of wave- vectors for all coupled waves. The ex- 
traordinary, distributed-feedback properties of OPA and 
thepossibility of BWMOPO in NIMs were predicted in 
0, [H, E3] ■ Herein, we explore the effects of construc- 
tive and destructive nonlinear quantum interference pro- 
cesses in NLO centers embedded in NIMs on fully and 
quasi-resonant four-wave mixing of the backward waves 



that were not investigated in 17J. Resonant coupling 
strongly enhances the NLO response of the composite. 
It appears that such a nanostructured composite exhibits 
extraordinary multiple-resonance properties with respect 
to the thickness of the slab, the density of the embedded 
centers, and the frequencies and strengths of the control 
fields. It also appears, that the linear phase-mismatch 
introduced by the host material can be negated at near- 
resonant FWM coupling of backward and ordinary waves. 
The investigation of such properties is important for the 
optimization of direct, coherent energy transfer from the 
control fields to the counter-propagating negative-index 
signal and the positive-index idler. This opens opportu- 
nities for the compensation of optical losses in NIMs and 
the creation of unique NIM-based photonic microdevices. 



BASIC IDEA 

The basic scheme of resonant FWM of the BW in a 
NIM is as follows. A slab of NIM is doped by four-level 
nonlinear centers [Fig. [lja)] so that the signal frequency, 
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FIG. 1: Scheme of quantum controlled FWM interaction (a) 
and coupling geometry (b). Here, 014 is frequency of the signal, 
uji and uz of the control fields, L02 of the idler, n(u4) < 0. 



0J4, falls in the NI domain \n(un) < 0], whereas all the 
other frequencies, oj\, u) 3 and lu 2i arc in the the positive 
index (PI) domain. Below, we show the possibility to 
produce transparency and even amplification above the 
oscillation threshold at 1V4 controlled by two lasers at wi 
and cj 3 [Fig. Hfb)]. These fields generate an idler at 
u>2 = CJ3 + lu\ — CJ4, which experiences population inver- 
sion or Raman amplification caused by the control fields. 
This opens an additional channel of energy transfer from 
the control fields to the signal. The amplified idler con- 
tributes back to u>4 — u>3+uji~ui2 through FWM and thus 
causes strongly enhanced OPA of the signal. Unlike ordi- 
nary off-resonant NLO, a many-order resonance enhance- 
ment of the NLO coupling is followed by a laser-induced 
strong change of the local optical parameters. The con- 
trol fields may cause population transfer and even inver- 
sion, as well as modulation of the probability amplitude, 
and split and quantum interference modifications of the 
resonance shapes. Through these effects, therefore, the 
control fields may serve as a tool for harnessing the local 
optical coefficients. Alternatively, such changes can be 
minimized so that the major amplification would come 
directly from the energy exchange between the control 
fields, the signal and the idler through the FWM pro- 
cesses. In fact, the interference of quantum pathways in 
the vicinity of the resonances may even lead to possibility 
that the overall process ceases to be seen as a set of suc- 
cessive one- and multi-photon elementary processes [HI 
(and references therein). 

First we shall show that the features of OPA in NIMs, 
which are employed in this work, appear in stark con- 
trast with the basic properties of such processes gener- 
ally known for ordinary PIMs. We assume that only uj 4 is 
close to the NI resonance. All other fields are far enough 
away to remain in the PI frequency domain. All wave 
vectors are presumed to be co-directed in order to en- 
sure maximum phase matching. Because the energy flow 
S4 = (c/An) [E x H] does not depend on the sign of the re- 
fractive index, it appears contradirected to k4 and, there- 



fore to all other wave vectors [Fig. UJ(b)]. This imposes 
extraordinary features on the nonlinear propagation of 
the coupled waves. The electric components of the waves 
and the nonlinear polarization are taken to be of the form 

Ej(z,t) = (E j /2)exp[i(k j z-bJ j t)]+c.c, (1) 
P* 2 L (z,t) = (P^ 2 L /2)exp{i[k 4>2 Z-u)^ 2 t]} + c.c. (2) 



Here, j = {1-4}, w 4 + w 2 = k^ 2 = ki + k 3 -k 2A , 

kj = KK/c > 0, P^ 2 L /2 = XilE! E 3 EZ a and X g 
are effective nonlinear susceptibilities. Then the slowly- 
varying amplitudes of the coupled waves at UI4 and w 2 
are given by the equations 



dEi/dz = -i^ v E^e^y[iAkz] + (a A /2)E A , 
dE 2 /dz = ij^ E'l exp[iAkz] - (a 2 /2)E 2 . 



(3) 
(4) 



Here, 7^4 = {Air\^ 2y A\ijjl A / k 2A c 2 )xi\EiE 3 are NLO cou- 
pling coefficients; Cj and fij are the dielectric permittiv- 
ities and magnetic permeabilities (which are negative at 
u>i); Ak = k\ + k 3 — k 2 — k A ; and ctj are the absorption or 
amplification coefficients. The amplitudes of the funda- 
mental (control) waves E\ and E 3 are assumed constant 
along the slab. Note that there are three fundamental dif- 
ferences in equations ([3] ) and ([4]) as compared with their 
counterparts in ordinary materials with positive perme- 
ability \l and permittivity e at all frequencies. First, the 

(3) 

sign of the nonlinear polarization term 74 is opposite 

(3) 

to that of 72 , which occurs because ^4 < 0. Second, 
the opposite sign appears for a A . Third, the boundary 
conditions for a A must be defined at the opposite side 
of the sample as compared to those for all other waves. 
These differences lead to a counterintuitive evolution of 
the signal and idler along the medium. Assuming £4.3 
constant and taking into account the boundary condi- 
tions E^{L) = E^l and E 2 (0) = (where L is the slab 
thickness), the solutions to equations (J3j- ((4j) are found 



Ei{z) = Aiexp[(#z]+A2exp[(#z], (5) 
E 2 (z) = K1A1 exp[/3j"z] + k 2 A 2 exp[/3^z], (6) 



where, 
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0i,2 ± (tAfc/2), /3i, 2 = (a 4 - a 2 )/(4) ± iR, 
M = E 4L k 2 /D, A 2 = -EiLKi/D, 
D = k 2 exp[/3i L] — ki exp\f3 2 L], 
= (±i? + is)/ 7 4, R 
7 2 74, s 



Kl,2 
„2 _ 



(a4 + a 2 )/4-t(A*/2). 



The transmission (amplification) factor for the negative- 



index signal, T±(z) = 
tion 

T 4 (0) = T 40 = 



\E A {z) / E 4L \ , is given by the equa- 



exp{- [(a 4 /2) -s]L} 



cos RL + (s / R) sin RL 



(7) 
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The fundamental difference between the spatial depen- 
dence of OPA in ordinary and NIM materials is explic- 
itly seen at <x, = Afc = in the off-resonant case, where 
72 = 74 = 9 an d hence g 2 is a real, positive quantity. 
Then equations ©-([7]) reduce to 

T 4 (z) = |cos( 5 z)/cos(.gL)| 2 , (8) 

whereas in ordinary media the signal would exponentially 
grow as T oc exp(2gz) Equation © and output 

magnitude T40 = |cos(<7-L)| 2 present a sequence of ge- 
ometrical resonances, with behavior similar to that of 
distributed-feedback resonance. Thus, in contrast with 
ordinary media, the spatial distribution of the signal in- 
side the slab and its output value at z = exhibit un- 
usual oscillatory dependence on the slab thickness. Such 
extraordinary resonances provide for the feasibility of 
attaining the oscillation threshold for the generation of 
the entangled counter-propagating left-handed, frw 4 , and 
right-handed, hw2, photons without a cavity. Such gener- 
ation would occur when gL — > (2j + 1)tt/2. The rate of 
the energy transfer from the control fields to the signal 
depends on the local strengths of both the signal and the 
idler. Thus, a change in their absorption and amplifica- 
tion rates my result in a dramatic change in distribution 
of the signal and the idler across the slab and, conse- 
quently, in their output values. Ultimately, the tailoring 
of the output signal stems from its extraordinary evolu- 
tion along the slab. 




(a) (b) 

FIG. 2: The effect of amplification of the idler and the dif- 
ference between the distribution of the signal, T 4 {z), and the 
idler, V2 = \E2{z) / 'E 4 {L)\ 2 , along the slab NIM, main plots 
(a); and in ordinary materials, (b), at otherwise identical con- 
ditions. Inset (a): geometrical resonances in output signal 
T 4 (z = 0). 1: a 2 L = a 4 L = 1, gL=1.9; 2: a 2 L = -3, 
a 4 L = 5, gL=11.049. 

An example of such dependence and its difference from 
that in the ordinary materials is depicted in Fig. [2a), (b). 
The inset in Fig. [2] (a) shows narrow resonances in the 
output signal at z — as the function of the quantity 
gL. The position and shape of the resonances depend 



on the absorption or amplification indices. The main 
plot shows the oscillation behavior of the signal and the 
idler across the slab at a slightly off-resonant values of 
gL, which results in the overall depletion of the signal 
and the idler. Since the idler grows toward the back 
facet of the slab at z = L and the signal experiences 
absorption or amplification in the opposite direction, the 
maximum of the signal at given parameters may appear 
closer to the back facet of the slab. As seen from the 
comparison of plots 1 and 2 in Fig.[2](a), such distribution 
depends strongly on the difference in the absorption or 
amplification rates for the signal and the idler. Hence, 
the changes in the slab thickness, or in the intensity of 
the control fields, or in their absorption indices lead to 
significant, resonance, changes in the output signal. Such 
a dependence is in stark contrast with its counterpart in 
positive-index materials shown in Fig. [2] (b) . 

RESONANT COUPLING AND QUANTUM 
CONTROL: NUMERICAL SIMULATIONS 

This main section is aimed at numerical experiments 
to show specific the counterintuitive features and data 
which are characteristic for resonant and near-resonant 
four- wave coupling of backward waves. The information 
presented here can be used for the all-optical tailoring of 
such coupling. 

Quasi-resonant coupling and compensation of phase 
mismatch 

An essential difference of the resonant nonlinear pro- 
cesses under investigation from ordinary off-resonant cou- 
pling is that the nonlinear susceptibilities 74 and 72 be- 
come complex and different from each other in this case. 
Hence, the factor g 2 may become negative or complex, 
which causes additional radical changes in the nonlinear 
propagation features. Figure [3] shows that resonant cou- 
pling and, hence, complex or negative magnitudes of g 2 
may lead to the result that the requirement Afc + 8k = 
ceases to be optimum for the process under investigation, 
which is also in striking contrast with the off-resonant 
OPA of the weak signal in ordinary materials. Some mis- 
match 5k, introduced by the host material, may bring 
the optimum signal frequency to the maximum of the 
NLO conversion rate. Alternatively, transparency and 
oscillation may become impossible for the frequencies 
that correspond to maximum nonlinear response unless 
the phase mismatch is elaborately added. Specifically, 
Fig. [3fa) and Fig. [3^b) prove the possibility of negating 
of the combine phase mismatch (Afc + 8k)L = (3/4)7r 
and (Ak + 8k)L = tt. Larger absolute values of negative 
Im(gL) 2 provide for compensation of larger positive val- 
ues of (Afc + Sk)L. Note, however, that the requirement 
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of the control fields are introduced as G\ = Eidi g /2h, 
G3 = E 3 d nm /2h. For the optical transitions indicated 
above, the magnitudes G ~ 10 12 s -1 correspond to the 
control field intensities in the range of / ~ 10 - 100 
kW/(0.1mm) 2 . The resonance offset is O4 = uj^ — ui m l> 
other resonance detunings Qj are defined in a similar way. 
The values X4,2 denote the effective linear susceptibilities, 
while «4o and Xio are their fully resonant values with all 
driving fields turned off. We assumed that 90% of the 
signal is absorbed by the host slab. 



FIG. 3: Correlation between the transmittance T4 of the NIM 
slab, Ke(gL), lm(gL) and combined linear phase mismatch 
introduced by the host material, 8k, and doping material, 
Afc. a 2 L = -3, a 4 L = 5. a: (Afc + 5k)L = (3/4)tt. (b): 
T4 = 1; solid line - phase mismatch (Afc + 8k)L — (3/4)-7r, 
dashed line - phase mismatch (Afc + 8k)L — n. 



of linear phase mismatch is known for coupling strong 
fields due to the fact that the ratio of the intensities of 
the coupled fields and, consequently, the coupling phase 
varies along the medium. A similar situation occurs here 
for the weak counter-propagating waves and the control 
fields that are homogeneous along the slab. 



Numerical model for embedded resonant NLO 
centers 

To demonstrate the outlined extraordinary NLO fea- 
tures and to prove the proposed possibilities of tailored 
transparency and oscillations, we have adopted the fol- 
lowing model for numerical simulations, which is char- 
acteristic for molecules imbedded in a solid host: the 
energy level relaxation rates are r„ = 20 x 10 6 s , 
120 x 10 6 



F — F — 

1 g x ra 



' 6 s 1 ; the partial transition prob- 



abilities are 7^ = 7 x 10° s , 7 gn 
7mn = 5 x 10 6 s-\ lml = 10 x 10 6 s- 1 
mogeneous transition half-widths are Ti g 
T lm = 1.9 x 10 12 s-\ F 
1.8 x 10 



■ ng 



= 4 x 10 b s-\ 
and the ho- 



r — 

nm — 
!l0 



^ „ = 5 x 10 io s- 1 ; T ln = 10 1 " s- 
We assume that A2 = 756 nm and A4 = 480 nm. The 
density-matrix method described in [18j | is used for cal- 
culation of the intensity-dependent local parameters, and 
we also included quantum nonlinear interference effects. 
This allows us to account for changes in absorption, am- 
plification and refractive indices as well as accounting 
for the striking changes in the magnitudes and signs of 
the NLO susceptibilities. The parameters Re 74,2, ^74^ 
and, hence, Re g 2 and Im g 2 may experience significant 
changes also due to the population redistribution over 
the coupled levels caused by the control fields. The lat- 
ter strongly depends on the ratio of the partial transi- 
tion probabilities. In the numerical experiments below, 
the coupling Rabi frequencies that represent the strength 



Fully resonant control fields 

Figure 2] depicts the local parameters entering in 
Eq. ([7]) with fully resonant control fields at the given 
strengths. All parameters are scaled to the indicated 
fully resonant values taken at all strong fields turned 
off. The quantity Rex4 represents changes in the re- 
fractive index at W4. Numerical analysis is done with 
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FIG. 4: Nonlinear interference structures in the local opti- 
cal parameters attributed to the embedded centers, y = 
Q 4 /(0.0ir !ro ). fij = fl 3 = 0, Gi = 1.646 GHz, G 3 = 
45.539 GHz. 

the aid of the steady-state solutions to the density ma- 
trix equations [18| . Here, y s = f24/(0.01F; m ). Modu- 



lation of the quantum amplitudes by the strong driv- 
ing fields and constructive and destructive interference 
of multistep and multiphoton quantum pathways at the 
coupled transitions cause the appearance of nonlinear 
interference spectral structures in the effective local lin- 
ear and nonlinear optical parameters [l8| (and references 
therein) . Figure [4] shows that the changes caused by the 
control fields in the local optical parameters may be com- 
mensurable with their original values. Consequently, the 
nonlinear propagation of the signal and idler can be tai- 
lored through such processes. Control fields at given in- 
tensities cause appreciable redistribution of the energy- 
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level populations accompanied by population inversion 
at the resonant idler transition: r; w 0.60, r g w 0.38, 
r n 0.01, r m m 0.01. Other optical parameters at 
LO4 = uj m i are as follow: a^/a^ ~ 0.11, a2/a 4 o « —0.18, 
Ak/a i0 ~ 0.28, Im(7 2 /a 40 ) » -0.008, Re(7 2 /a 40 ) = 0, 
Im( 74 /a4o) w -0.002, Re( 74 /a 40 ) = 0, Im( 5 2 /al ) = 0, 

Re(g 2 /a 40 ) ~ 1-6 x 10~ 5 . 

Figure [5] displays the resonance dependence of the out- 
put signal on the optical thickness of the slab, on the res- 
onance offset of the signal, on the strengths of the control 
fields, and the distribution of the signal and idler along 
the slab as tailored by the given control fields. As dis- 
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FIG. 5: Geometrical resonances in the output signal and idler, 
(a) and (b); dependence of the output signal on the intensity 
of the control fields at L=L r , (c); and distribution of the signal 
and the idler inside the slab at L slightly different from L r , 
(d). y = n 4 /(0.0ir !m ). Qj = Q 3 = 0, a NIM L = 2.3. (a), (c) 
and (d): uj 4 = w im . (a), (b) and (d): Gi = 1.646 GHz, G 3 = 
45.539 GHz. (c): L r /L ra = 150.3130, single peak (dash- 
dot) - Gi = 1.655 GHz, solid line - Gi = 1.657 GHz, dashed 
line - Gi = 1.678 GHz; (d): L/L ra = 150.31. 

cussed above, the output signal for the waves coupled in 
the NIM slab through OPA exhibits a set of distributed 
feedback-type resonances. Such resonant behavior ap- 
pears with respect to the intensity of the fundamental 
fields, as well as with respect to the product of the slab 
length and the density of the nonlinear centers (the op- 
tical density), and also with respect to the resonance off- 
sets for the signal and fundamental fields. Optimization 
of the output signal at z — is determined by the in- 
terplay of absorption, idler gain, FWM and, hence, by 
the wave vector mismatch. This appears to be a multi- 
parameter problem associated with sharp resonance de- 
pendencies. Figure O displays the results of the numerical 
analysis of Eq. ([7]) based on the steady-state solutions to 
the density-matrix equations for one of the transmission 



resonances. In the given case, the strongest and most 
easily achieved resonance corresponds to the center of 
the signal transition at w 4 = w TO ; . Here, we introduce the 
scaled product of the slab length and the number density 
of the embedded centers, L/L ra , through the resonance 
absorption length, L ra = . Figure [5]^a) displays a 
narrow geometrical resonance as a function of the slab 
thickness or the density of the embedded NLO centers 
at L = L r = 150.3130 L ra for the optimum frequency 
offset for the signal at Q4 — 0. Figure [H^b) shows that 
the width of the transparency window is on a scale less 
than the narrowest (here Raman) transition half-width 
and the resonant absorption length. Figure [5f c) shows 
that the output signal at the given frequency and optical 
density of the slab is strongly dependent on the magni- 
tude and ratio of the control field intensities. Figure [5jd) 
shows that the signal intensity inside the slab may sig- 
nificantly exceed its output value at z=0, and the idler 
intensity substantially exceeds that for the signal. Such 
features depend on the ratio of the OPA, absorption rate 
of the signal and amplification rate of the idler. Here, 
the slab thickness (or density of the nonlinear centers) 
is taken to be slightly different from its resonance value. 
Calculated 



Frequency offsets on the order of the narrowest 
Raman transition width 
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FIG. 6: Nonlinear interference structures in absorption of the 
signal and amplification of the idler, (a); in refractive indices 
and phase mismatch attributed to the embedded centers, (b) 
and (c); and in the nonlinear coupling parameters, (d)-(f). 

y s = n 4 /r ln . fii = -n 3 = r in , d = 21 gh z , g 3 = 27 gh z . 



In this case, nonlinear interference structures with dif- 
ferent spectral widths only partially overlap each other. 
Typical results for the numerical experiments are given 
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below. Population transfer, power broadening of the 
nonlinear interference structures, and the overall power- 
saturation of the local optical and NLO parameters de- 
pend on both the intensities and resonance frequency 
offsets of the control fields. This determines different 
sets for the optimum parameters for coherent energy ex- 
change between the control fields, backward signal and 
the idler. Figure [6] displays an appreciable change in the 
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FIG. 7: Geometrical and frequency resonances in the output 
signal and idler, (a) and (b); and distribution of the signal 
and idler inside the slab at L slightly different from L r , (c). 

y B = n 4 /r ln . fii = -n 3 = r !n , d = 21 gh z , g 3 = 27 GHz, 

Q-nimL = 2.3, 8k/ct4o = —0.013069. (a) and (c): y S Q — 
-5.624. (b): L r /L ra = 41.86 (c): L/L ra = 41.865 

local optical parameters at the given frequency detunings 
and the optimum intensity of the control fields as com- 
pared with Fig. [U Consequently, the propagation fea- 
tures change such that the maximum transparency and 
amplification of the signal become achievable at a slab op- 
tical density about three times smaller L r = 41.86 a^Q 1 
[Fig. E[a)]. Other features are similar to those of Fig. 
HI Energy-level populations at given strength and reso- 
nance offsets of the control fields and other local optical 
parameters at y s = y s o are as follow: r; « 0.49, r g m 0.48, 
r n ph 0.015, r m rs 0.015, a 4 /a 4 o ~ 0.38, a 2 /a 4 o ~ —0.36, 
Afc/a 40 w 0.15, lm(7 2 /o;4o) « -0.046, Re(7 2 /a 4 o) ~ 
-0.037, Im(7 4 /a 40 ) w 0.059, Re(7 4 /a 40 ) « -0.047, 
lm(g 2 /al ) rj -4.26 x lO" 3 , Re(g 2 /a% ) ~ ~ 8 - 4 x 10 ~ 4 - 



Frequency offsets on the order of the optical 
transition width 



Here, optimum coupling parameters are analyzed for 
frequency resonance offsets that are several hundred 
times larger than those in the previous subcase. Figure[8] 
displays the appearance of two nonlinear resonances and 
an increase in the amplification of the idler. At the given 
parameters, the control fields cause essential change in 
the level populations, n ~ 0.49; r g w 0.48; r n « 0.015; 
r m 0.015, which is accompanied by the large popu- 
lation inversion at the idle transition. Figure [9] shows 
the transmitted signal at the frequency close to the max- 
imum NLO response of the embedded molecules, which 
may become possible owing to the appropriate phase mis- 
match, Sk, introduced by the host material. Figure a) 
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FIG. 8: Scaled absorption index for the signal and amplifica- 
tion index for the idler (a); changes in the refraction indices 
(b) and phase mismatch (c); NLO coupling parameters, (d) 
(e); and parameter g 2 , (f), in the presence of the control fields 
at Gi = 33.19 GHz, G 3 = 15 GHz, Qi = -3r 5 Z, Q 3 = 0. 
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FIG. 9: Geometrical and frequency resonances in the output 
signal and idler, (a) and (b); and distribution of the signal 
and the idler inside the slab at L slightly different from L r , 
(c). j/4 = n 4 /r mi , Qi = -3Fgl, Q 3 = 0, Gi = 33.19 GHz, 
G 3 = 15 GHz. a NIM L = 2.3, y 40 = -1.578947, 6k/a 40 = 
3.165 x 10" 4 . (a) and (c): y 4 = y w . (b): L r /L ra = 40.9977. 
(c): L/Lra =40.8977. 



displays a narrow geometrical resonance at L w 41L ra 
for the optimum frequency offset for the signal at f2 4 =- 
1.578947 r m /. This provides the optimum compensation 
the given phase mismatch, 8k, introduced by the host 
material. Changes in Afc introduced by molecules are 
accounted for within the simulations. Figures [9][b) and 
(c) depict qualitatively similar behavior as in the previous 
subcases. The calculated optical parameters for ?y 4 = y 4 o 
are as follow: Im(7 2 /a 4 o) ~ —0.014, Re^/a^o) ~ 
-0.015, Im(7 4 /a 40 ) « 0.017, Re( 74 /a 40 ) ~ -0.019, 
Im( 5 2 /al ) « -5.31 x 10~ 4 , Re( 5 2 /ai ) = 6.25 x 10~ 5 , 
a 4 /a 4 o = 0.15, a 2 /a 4 o = —0.42, A/c/a 4 o = —0.17. 

Amplification in the output signal maxima in Figs. [H 
[6] and [8] reaches many orders of magnitude, which indi- 
cates the feasibility of oscillations and, hence, cavity-less 
generation of counter-propagating signal and idler pho- 
tons. It is known that even weak amplification per unit 
length may lead to lasing provided that the correspond- 
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ing frequency is in a high-quality cavity or feedback res- 
onances, which is equivalent to a great extension of the 
effective length of a low-amplifying medium. Assuming 
a resonance absorption cross-section 040 ~ 10 -16 cm 2 , 
which is typical for dye molecules, and a concentration of 
molecules N ~ 10 19 cm' 3 , we obtain ~ 10 3 cm' 1 , 
and the required slab thickness is in the range L ~ (10 — 
100) /a. The contribution to the index of refraction by the 
impurities is estimated as An < 0.5(A/47r)a4 ~ 10~ 3 , 
which essentially does not change the negative refractive 
index. 



CONCLUSION 

In conclusion, we have investigated the feasibility of 
all-optical manipulation of the optical properties of NIMs 
through coherent nonlinear-optical energy transfer from 
the control to the signal field. The strong nonlinear op- 
tical response of the composite is primarily determined 
by the embedded resonant four-level nonlinear centers 
and, hence, can be adjusted independently. In addi- 
tion, we have shown the opportunity for quantum con- 
trol of the local optical parameters, which employs con- 
structive and destructive quantum interference tailored 
by two auxiliary driving control fields. Such a possibil- 
ity is proven with the aid of a realistic numerical model. 
The investigated features are promising for the compen- 
sation of losses in strongly absorbing NIMs, which is 
the key problem that limits numerous revolutionary ap- 
plications of this novel class of electromagnetic materi- 
als. Among the other possible applications are a novel 
class of the miniature frequency-tunable narrow-band fil- 
ters, quantum switchers, amplifiers and cavity-free micro- 
scopic optical parametric oscillators that allow the gen- 
eration of entangled counter-propagating left- and right- 
handed photons. The unique features of the proposed 
photonic devices are revealed, such as the strongly reso- 
nant behavior with respect to the material thickness, the 
density of the embedded resonant centers and the inten- 
sities of the control fields, the feasibility of negating the 
linear phase-mismatch introduced by the host material, 
and the important role of the supplementary nonpara- 
metric amplification of the idler. 
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